Limited phylogeographic and genetic connectivity in Acacia species of low stature in an arid landscape

Abstract Widespread plant species are expected to maintain genetic diversity and gene flow via pollen and seed dispersal. Stature is a key life history trait that affects seed and potentially pollen dispersal, with limited stature associated with limited dispersal and greater genetic differentiation. We sampled Hill’s tabletop wattle (Acacia hilliana) and curry wattle (Acacia spondylophylla), two co‐distributed, widespread, Acacia shrubs of low stature, across the arid Pilbara region of north‐western Australia. Using chloroplast sequence and nuclear microsatellite data we evaluated patterns of population genetic and phylogeographic diversity and structure, demographic signals, ratios of pollen to seed dispersal, evidence for historical refugia, and association between elevation and diversity. Results showed strong phylogeographic (chloroplast, G ST = 0.831 and 0.898 for A. hilliana and A. spondylophylla, respectively) and contemporary (nuclear, F ST = 0.260 and 0.349 for A. hilliana and A. spondylophylla, respectively) genetic structure in both species. This indicates limited genetic connectivity via seed and pollen dispersal associated with Acacia species of small stature compared to taller tree and shrub acacias across the Pilbara bioregion. This effect of stature on genetic structure is superimposed on moderate levels of genetic diversity that were expected based on widespread ranges (haplotype diversity h = 25 and 12; nuclear diversity He = 0.60 and 0.47 for A. hilliana and A. spondylophylla, respectively). Contemporary genetic structure was congruent at the greater landscape scale, especially in terms of strong genetic differentiation among geographically disjunct populations in less elevated areas. Measures of diversity and connectivity were associated with traits of greater geographic population proximity, population density, population size, and greater individual longevity, and some evidence for range expansion in A. hilliana. Results illustrate that low stature is associated with limited dispersal and greater patterns of genetic differentiation for congenerics in a common landscape and highlight the complex influence of taxon‐specific life history and ecological traits to seed and pollen dispersal.


| INTRODUC TI ON
Understanding drivers of species genetic diversity, differentiation and demography is of central interest to the fields of biogeography, phylogeography, population and conservation genetics, and ecology (Avise, 2000). Species with sympatric regional biogeographic distributions are expected to have experienced the same historical geographic, geological, and/or climatological phenomena across landscapes (Arbogast & Kenagy, 2001). As a result, comparative phylogeographic analyses of broadly co-distributed species can assist in elucidating the impact of regional geological and other landscape features, as well as historical macroevolutionary climatic and environmental events, on patterns of genetic diversity, differentiation, and demography at the landscape scale (Allendorf et al., 2013;Arbogast & Kenagy, 2001;Avise, 2000). Identifying common patterns in areas of haplotype diversity and divergence can lead to the identification of landscape features that may have provided potential historical refugia, or barriers to dispersal, as well as patterns of population persistence, contraction, and/or expansion (Hickerson et al., 2010;Keppel et al., 2012).
Identifying macroevolutionary patterns of biogeographic concordance among co-distributed taxa should also be balanced by biologically informed approaches that determine the contribution of taxon specific life history and ecological traits on levels and structuring of genetic diversity on more microevolutionary temporal and spatial scales (Papadopoulou & Knowles, 2016;Sullivan et al., 2019).
For example, plant species are influenced to varying degrees by a range of life-history traits that drive fine-scale processes of seed and pollen dispersal, as well as by ecological traits (Broadhurst et al., 2017;Hamrick & Godt, 1996). Key life history traits that affect seed and pollen dispersal include mating systems (Duminil et al., 2007), the mechanisms of dispersal (Duminil et al., 2007;Thomson et al., 2010), and stature (Petit & Hampe, 2006;Thomson et al., 2010). Meta-analyses have shown these key traits, as well as geographic and demographic attributes of range size, degree of population geographic disjunction, localized abundance or density, population size, and individual longevity and fecundity, have the greatest effects on genetic diversity and differentiation (Broadhurst et al., 2017;Gitzendanner & Soltis, 2000;Hamrick, 1979;Hamrick & Godt, 1996;Hamrick et al., 1992;Loveless & Hamrick, 1984;Nybom, 2004;Nybom & Bartich, 2000). The use of comparative phylogeography in deducing how differing life history and ecological traits contribute to macro and microevolutionary patterns of genetic diversity, differentiation, and demography of specific plant taxa is not common (although see García-Verdugo et al., 2019;Massatti & Knowles, 2014;Paz et al., 2015). Comparative approaches that assess concordance in contemporary patterns of genetic diversity driven by trait based microevolutionary processes, as well as historical abiotic macroevolutionary processes, offer fundamental insights into the taxa themselves and allow for refinement of predictions on how taxon-specific traits affect diversity, differentiation, and demography (Papadopoulou & Knowles, 2016).
Historical evolution in a common landscape is typically expected to result in concordant biogeographic patterns (Arbogast & Kenagy, 2001;Hewitt, 2004). However, the influence of intrinsic taxon-specific life history and ecological traits on phylogeographic patterns, as well as contemporary population genetic diversity and differentiation, is often underestimated (Papadopoulou & Knowles, 2016). Biogeographic studies are also underrepresented in many species-rich environments, particularly in the Southern Hemisphere (Beheragaray et al., 2008;Hickerson et al., 2010), in Australia (Beheragaray et al., 2008), and in the world's arid zones Lexer et al., 2013;Robertson & Chan, 2011). The Pilbara bioregion of north-western Australia provides an opportunity to investigate comparative genetic patterns in a species-rich, topographically complex, arid environment in the Southern Hemisphere.
Previous studies of three Acacia present in the Pilbara have shown evidence for congruent patterns of historical population isolation and persistence at the broad regional scale (Levy et al., 2016;Nistelberger et al., 2020). However, patterns of phylogeographic and contemporary population genetic diversity and differentiation have also been shown to vary with differences in taxon-specific life history and ecological traits known to impact pollen and seed dispersal, including differences in species geographic range size, and growth form or stature.
For example, species with widespread ranges, A. ancistrocarpa and A. pruinocarpa, show greater levels of chloroplast haplotype diversity, haplotype network complexity, phylogeographic structure, and contemporary genetic diversity and connectivity (Levy et al., 2016;Nistelberger et al., 2020), compared to the range restricted A. atkinsiana (Levy et al., 2016). Also, the tallest of the three species, A. pruinocarpa (to 12 m), has increased haplotype diversity, a more complex haplotype network, and some evidence for population expansion (Nistelberger et al., 2020) in comparison to the smaller shrub species (both to 4 m) that show more limited haplotype network diversity, and evidence for population persistence (Levy et al., 2016).
Limiting comparisons to congeneric species that share key life history traits that primarily affect dispersal allows for refined hypotheses regarding the contributions of life history and ecological traits to diversity, differentiation, and demography. Here, we test hypotheses based on intrinsic taxon-specific life history and ecological traits to improve understanding of macro and microevolutionary drivers of genetic diversity and differentiation for Acacia species in the Pilbara. We use both chloroplast (cpDNA) sequencing data and nuclear microsatellite (nDNA) data to assess intraspecific patterns of phylogeographic and population genetic diversity and structure in Hill's tabletop wattle (Acacia hilliana Maiden, Figure 1a) and curry wattle (Acacia spondylophylla F. Muell. Figure 1b). Acacia hilliana and K E Y W O R D S Acacia, arid zone, genetic connectivity, genetic diversity, Pilbara, stature A. spondylophylla are regionally common, with widespread, broadly co-distributed geographic ranges within the Pilbara (Figure 1). Both species have widespread ranges extending outside the Pilbara region suggesting that haplotype and nuclear diversity will be comparable to similarly widespread Acacia species (Levy et al., 2016;Nistelberger et al., 2020). However, unlike the Acacia species studied so far, which are large shrubs or trees, A. hilliana and A. spondylophylla are both small, flat-topped shrubs less than 1.5 m tall. A global review has shown that low stature is associated with significant limitations to seed dispersal distance for species with the same seed dispersal syndrome (Thomson et al., 2011). Given that Acacia in the arid Pilbara may be expected to have similar key life history traits that primarily affect pollen and seed dispersal, a predominantly outcrossed mating system, pollination by generalist insects and seed dispersal primarily by gravity and ants, low stature may be expected to result in greater limitation to seed dispersal and greater phylogeographic structure for the study species. The two study species differ in other life history and ecological traits including the degree of population geographic disjunction, localized population density, population size, and individual longevity and fecundity (detailed in the Materials and Methods), that may also result in observable differences in contemporary patterns of genetic diversity and differentiation between the species.
We use a standardized comparative approach at macro and microevolutionary spatial and temporal scales to gain insight into how levels and patterns of genetic diversity are influenced by taxonspecific life history and ecological traits. Our assessment of two widespread Acacia shrubs of small stature is placed in the context of highly comparable studies, that used a similar suite of chloroplast sequencing data and nuclear microsatellite data and the same sampling approach and sampling numbers, on several Acacia with variable species ranges across the Pilbara and that are either trees or large shrubs. We make the following predictions for A. hilliana and A. spondylophylla; (i) they will exhibit levels of haplotype and nuclear diversity and contemporary genetic connectivity comparable to those of other widespread Acacia present in the Pilbara; (ii) given that low growing habit likely leads to more restricted seed dispersal, they will exhibit stronger phylogeographic structure compared to that of tall F I G U R E 1 Geographic distribution of sampled populations, chloroplast DNA haplotypes, and evolutionary relationships amongst chloroplast DNA haplotypes of (a) Acacia hilliana and (b) Acacia spondylophylla in the Pilbara bioregion. Populations are represented by pie charts. Segment size within each pie-chart corresponds to the proportion of individuals with that haplotype. Population codes correspond to those in Table 1. Greyscale indicates elevation. Top inset shows species known geographic distribution as blue circles. Pilbara bioregion indicated. Bottom inset Median-Joining networks show the evolutionary relationships between cpDNA haplotypes in sampled populations. Haplotypes are represented by colored circles. Color of circle corresponds to haplotypes in the geographic map. Circle size is proportional to the number of individuals with that haplotype and branch lengths are proportional to the number of mutations with lengths greater than one indicated by cross bars. Small black circles indicate unsampled hypothetical haplotypes. Panels to the right show the prostrate nature of representative individuals shrub and tree Acacia present in the Pilbara. Patterns of contemporary genetic structure and diversity are likely to be largely congruent for A. hilliana and A. spondylophylla as they have experienced the same historical geographic, geological, and climatological phenomena at the regional scale, although a greater degree of population disjunction, lower localized population density, smaller populations, and shorter individual longevity and lower fecundity in A. spondylophylla leads to another prediction of (iii) decreased genetic diversity within populations and increased genetic variation among populations of A. spondylophylla. We also assessed phylogeographic concordance for A. hilliana and A. spondylophylla in terms of the major landscape features of the Hamersley and Chichester Ranges, that may act as potential historical refugia, and the impact of population elevation on measures of contemporary diversity.

| Study region
The Pilbara bioregion of Western Australia is a highly structured landscape, recognized as a biogeographic transition zone between the tropics to the north, the sandy deserts to the east and the arid rangelands to the south (Cracraft, 1991;Pepper et al., 2013). The region supports a wide range of habitats with a climate dominated by high temperatures and limited annual rainfall, although localized flooding can result from thunderstorms that often follow tropical summer cyclones (Leighton, 2004). and tableland with jagged peaks, gorges, rolling hills, and tree-lined water courses with spinifex, hummock grasslands and Eucalyptus leucophloia steppe. The region is recognized as a significant center of species endemism across the flora and fauna biota of the Australian continent (Burbidge et al., 2010;Cracraft, 1991;Doughty et al., 2011;Durrant et al., 2010;Eberhard et al., 2005;Gibson & McKenzie, 2009;Guthrie et al., 2010;Ladiges et al., 2006;McKenzie & Bullen, 2009;Unmack, 2001). The Pilbara is also an important center of diversity for the Acacia genus (Maslin et al., 2010), whose species are common in Australian arid and semi-arid zone plant communities (Beadle, 1981).

| Study species
Hill's tabletop wattle (A. hilliana) is a semi prostrate sub-shrub reaching 0.3-1 m in the Pilbara. Plants have bright yellow, erect, cylindrical flower spikes to 45 mm from April to November (Maslin et al., 2010).
Brown seeds up to 5.5 mm form in erect, woody, resinous pods that produce a light but distinct citrus odor (Maslin et al., 2010). Each seed has a large eliasome and seed is shed from pods following dehiscence. Individuals are slow growing, expected to live 50+ years, and are highly fecund (seed production). The species has a widespread but scattered distribution from the Pilbara through central and northern Australia to north-western Queensland (Figure 1a).
Plants occur on skeletal red loamy sand or sand over laterite, quartzite, or ironstone, on sand dunes and on spinifex plains in deep sand on the lowlands, but also on rocky elevated sites on the ranges and plateaus of the Pilbara (Maslin et al., 2010). Populations are generally large, covering many hectares, and are continuous and locally dense on lowland rolling hills and upland plains of the ranges, although populations off the ranges are less continuous with some disjunctions (S van Leeuwen, personal observation).
Curry wattle (Acacia spondylophylla) grows slightly larger but is still a small, spindly shrub typically 1-1.5 m. Plants have resinous leaves, stems, and pods that produce a curry-like odor, with larger, bright yellow, round inflorescences appearing from May to August (Maslin et al., 2010). Slightly smaller black seeds, 3.5-4 mm long, with a small eliasome form within flat pods and are easily shed from pods (Maslin et al., 2010). Individuals are fast growing, living 5-15 years, and fecundity is limited with relatively few seeds being produced. The species has a scattered distribution across the Pilbara and central Australia (Figure 1b), with typically open and sparse populations that are more disjunct, less dense, and smaller than those of A. hilliana (S van Leeuwen, personal observation). Plants grow in shallow, sandy or rocky soil on ironstone hills and in gullies, as well as along creeks and colluvial flats on spinifex plains (Maslin et al., 2010).
Populations are patchily distributed or disjunct and occur with higher density on the Hamersley Range, although elevated populations are still more scattered, less dense, and smaller than those of A. hilliana (S van Leeuwen, personal observation).
These species are expected to have the same key life history traits as other Pilbara Acacia studied so far. Although there are exceptions, most Acacia share life history traits that have primary effects on dispersal of both seed and pollen. These include pollen dispersal by generalist insect pollinators (Bernhardt, 1989;Stone et al., 2003), self-incompatibility, or at least partial selfincompatibility that may be prezygotic or postzygotic (Kenrick & Knox, 1989;Morgan et al., 2002;New, 1984), and seed dispersed primarily by gravity and secondarily by myrmecochory (Berg, 1975;Davidson & Morton, 1984).

| Sampling and DNA extraction
Distributional information on each species across the Pilbara bi-
Variable sites were scored as a single mutation and indels coded as single binary states regardless of length. Indels arising from mononucleotide repeats were removed. These sequence data are available at the GenBank (www.ncbi.nlm.nih.gov/genbank) database under accession numbers MW138725-MW138756 and MW696776-MW696791.

| Nuclear microsatellite library construction, primer assessment, and locus assessment
Species specific nuclear microsatellite markers were developed for each species through partial genome sequencing on a 454 platform at the AGRF, using the methods of Gardner et al. (2011).
The resulting microsatellite libraries were sorted for functionality and the best 40 primers tested for polymorphism in a subset of individuals from multiple populations. Primer pairs were optimized for 16 reliable microsatellite loci for each species (Table S1).
Microsatellite amplification was performed in 7.25 μl multiplexed reactions containing 3.5 μl QIAGEN Multiplex PCR Master Mix (Qiagen), 0.75 μl primer mix, 2 μl sterile distilled water and 1 μl (10 ng) of template DNA. Cycling conditions were the same for all loci for both species; 95°C for 15 min, followed by 35 cycles of 94°C for 30 s, 60°C for 90 s and 72°C for 60 s, followed by a final extension of 60°C for 30 min. Forward primers contained a fluorescent color of the G5 label set (FAM, VIC, NED or PET). One μl of diluted PCR product was added to 12 μl of GeneScan™ LIZ®500 (-250) size Standard (Applied Biosystems)/formamide, and PCR products were visualized on a Biosystems 3730 Sequencer (Applied Biosystems). Genotypes for 20 individuals from each population were scored using Genemapper™ v.3.7 (Applied Biosystems) software with amplification and fragment analysis repeated once for any individuals that failed to amplify or produce scorable bands in the first instance. Allele bins were manually assigned and automatically checked. Sequences containing the assayed microsatellites are available at the GenBank (www.ncbi.nlm.nih.gov/genbank) database under accession numbers MTO80657-MTO80688.
Utility of loci was assessed for each species. Departure from Hardy-Weinberg equilibrium (HWE) and Linkage Disequilibrium (LD) among locus pairs were tested using exact tests implemented in GENEPOP v4.2 (Raymond & Rousset, 1995). Frequency of null alleles and their effect on differentiation estimates was assessed with FREENA (Chapuis & Estoup, 2007). Significant (p < .0001 with Bonferroni correction applied) departure from HWE was found in 15 out of 352 locus/population combinations for A. hilliana and in nine out 240 locus/population combinations for A. spondylophylla, but never for more than four loci per population. There was potential evidence of null alleles at some locus/population combinations for both species. Null allele frequencies ranged from 0.00 to 0.376 for A. hilliana and from 0.00 to 0.365 for A. spondylophylla (Table S2). However, F ST values were not significantly different in raw and corrected data sets (data not presented) so null alleles are likely to have little effect on related analyses (Chapuis & Estoup, 2007). Significant (p < .0001 with Bonferroni correction applied) LD was detected in 14 out of 120 locus combinations for A. hilliana and in two out of 120 locus combinations for A. spondylophylla.
Significance of values were assessed using a null distribution of 10,000 coalescent based simulations. Measures of phylogenetic diversity within populations (h S ) and total diversity (h T ), as well as coefficients of phylogenetic differentiation for unordered (G ST ) and ordered alleles (N ST ), were calculated with PERMUT v2.0 following Pons and Petit (1996). Phylogeographic structure was evaluated by assessing whether N ST was significantly greater than G ST using the U test for 1000 permutations in PERMUT.
spondylophylla was highly differentiated from all other populations due to a large number of unique alleles in high frequency. It is suspected that this population represents the rare hybrid Acacia adoxa var. adoxa × spondylophylla that occurs at low frequencies in the Ophthalmia Range west of Newman and from near Nullagine where the two parent species are common (Maslin et al., 2010). As a result, the Newman population that was collected as A. spondylophylla was removed from further analyses of the nuclear and chloroplast genomes.
Genetic structure was investigated using the Bayesian assignment approach in structure v2.3.2.1 (Pritchard et al., 2000). Analysis used the admixture ancestry model with the assumption of correlated allele frequencies amongst samples given the likelihood of high population similarity due to gene flow (Falush et al., 2003).
Known population locations were used as priors. A burn-in period of 10,000 was applied with 100,000 MCMC replications to assess K values ranging from K = 1 to K = the number of populations plus one, with 10 iterations of each K value. The optimal K value was de-

| Chloroplast haplotype diversity
For A. hilliana, 2368 cpDNA sites were scored. There were 33 variable sites with 30 parsimony-informative sites, leading to the identification of 25 haplotypes (

| Nuclear diversity and differentiation
Nuclear microsatellite diversity was moderate for both species (  In structure analysis of A. spondylophylla, L(K) increased steadily leveling off at K = 5 but with a large standard deviation, then continued to increase leveling off with no significant deviation at K = 11 TA B L E 2 Estimates of sequence diversity parameters, divergence parameters and neutrality tests for three and four chloroplast regions in Acacia hilliana and Acacia spondylophylla in populations from the Pilbara bioregion. Diversity and divergence parameters include total number of variable sites (S), the total number of haplotypes (h), overall haplotype diversity (Hd), nucleotide diversity ( Table 1. (d) Geographic distribution of genetic clusters. Colors correspond to those in (c). Sections of pie charts show the proportion of assignment of the given population to each genetic subcluster. Labels correspond to population codes in Table 1 Measures of genetic divergence among populations were high for each species. Genetic fixation among populations was lower for

| DISCUSS ION
Comparative analysis of congeneric species provides a phylogenetically controlled evaluation of how taxon-specific life history and ecological traits shape genetic diversity and connectivity.
Our analysis of two Acacia species showed levels of genetic diversity comparable to those of other co-distributed, widespread Acacia species in the arid Pilbara landscape, confirming expectations of the effect of species range on levels of genetic diversity.
Evaluation showing greater phylogeographic structure for the two low-statured study species, compared to other Pilbara Acacia species that are either large shrubs or trees, was consistent with the predicted association of stature with limited seed dispersal, as well F I G U R E 3 Assignment of individuals and populations to genetic clusters identified via Bayesian analysis of multilocus nuclear microsatellite genotype data for Acacia spondylophylla in the Pilbara bioregion. (a) Mean values of the ad hoc test statistic ln probability of the data (L(K)), values of the ad hoc test statistic DeltaK (ΔK), and bar plot indicating the assignment of individuals and populations to one of two optimally determined genetic clusters. Each individual is represented as a vertical line partitioned into two segments whose length is proportional to the individual coefficients of membership in each of two clusters. (b) Hierarchical clustering analysis of populations from Cluster One (Orange) in (a). Mean values of the ad hoc test statistic ln probability of the data (L(K)), values of the ad hoc test statistic DeltaK (ΔK), and bar plot indicating the assignment of individuals and populations to one of two optimally determined genetic clusters. Each individual is represented as a vertical line partitioned into two segments whose length is proportional to the individual coefficients of membership in each of two clusters. (c) Hierarchical clustering analysis of populations from Cluster Two (Blue) in (a). Mean values of the ad hoc test statistic ln probability of the data (L(K)), values of the ad hoc test statistic DeltaK (ΔK), and bar plot indicating the assignment of individuals and populations to one of nine optimally determined genetic clusters. Each individual is represented as a vertical line partitioned into nine segments whose length is proportional to the individual coefficients of membership in each of nine clusters. (d) Geographic distribution of eleven genetic clusters from the hierarchical analysis. Colors correspond to those in (b) and (c). Sections of pie charts show the proportion of assignment of the given population to each genetic subcluster. All population codes correspond to those in Table 1 as limited pollen dispersal. Also, as predicted, while patterns of genetic structure were largely concordant at the landscape scale, more subtle differences in diversity, differentiation and demography were evident and consistent with differences in the life history traits of these species.

| Patterns of diversity and phylogeographic structure
Both haplotype and nuclear diversity in A. hilliana and A. spondylophylla were comparable to those found for two other widespread F I G U R E 4 Principal co-ordinates analysis of the genetic distance between sampled individuals of (a) Acacia hilliana and (b) Acacia spondylophylla in the Pilbara bioregion. Population codes correspond to those in Table 1 F I G U R E 5 Heatmaps of pairwise population genetic divergence (F ST ) for nuclear microsatellite data for (a) Acacia hilliana and (b) Acacia spondylophylla in the Pilbara bioregion. Population codes correspond to those in Table 1 Acacia species with populations in the Pilbara bioregion, A. ancistrocarpa (Levy et al., 2016) and A. pruinocarpa (Nistelberger et al., 2020).
This result supports our first prediction regarding genetic diversity and is consistent with population genetic theory that plant species with widespread ranges maintain greater diversity than species with restricted range sizes, as they are expected to be less impacted overall by stochastic environmental and genetic effects that may negatively affect population sizes (Allendorf et al., 2013;Ellstrand & Ellam, 1993;Godt & Hamrick, 1993). The levels of haplotype and nuclear diversity in these widespread Acacia contrasts with lower haplotype and nuclear diversity detected for the range limited and near-endemic A. atkinsiana (Levy et al., 2016). Range size or distribution is known to have a strong effect on genetic diversity of the Australian flora (Broadhurst et al., 2017). These results confirm a significant influence of species geographic range on the maintenance of genetic diversity over both macro-and micro-evolutionary temporal scales in this arid landscape.
Detection of significant patterns of phylogeographic structure also matched our predictions. Seed dispersal is a key life history trait affecting gene flow among populations, with significant phylogeographic structure in chloroplast haplotypes indicating historical limitations to gene flow via seed dispersal. This would be expected for Acacia in the arid Pilbara, as primary dispersal of seeds is by gravity and secondary dispersal by myrmecochory. Gravity dispersal is likely to result in short distance primary seed dispersal of a few meters, while longer distance dispersal events are possible via ant dispersal (a few hundred meters, Pascov et al., 2015), although ant dispersed species have also shown some of the shortest dispersal distances documented, of just a few meters (Gómez & Espadaler, 1998;Hughes et al., 1994;Hughes & Westoby, 1992). For species that lack longer distance seed dispersal mechanisms, significant limitations to seed dispersal are most strongly associated with stature, with species of low stature having more restricted seed dispersal distances (Duminil et al., 2009;Nybom, 2004;Thomson et al., 2011). Our result of more pronounced phylogeographic structure in A. hilliana and A. spondylophylla than that found for populations of both the widespread large (to 4 m) shrub A. ancistrocarpa (G ST = 0.654, Levy et al., 2016) and the tall (to 12 m) tree species A. pruinocarpa (G ST = 0.230, Nistelberger et al., 2020), is consistent with this prediction and suggests that growth form or stature is a key life history trait driving seed dispersal distance in these Acacia. The lack of detectable phylogeographic structure found in A. atkinsiana (Levy et al., 2016), which is also a medium sized (to 3.5 m) shrub but with more restricted distribution, is more likely to be associated with the maintenance of a very limited number of extant haplotypes (four) over a shorter geographic range.
While we observed broad concordance in genetic patterns with our study species following predictions based on distribution and stature, we also observed differences between them in terms of contemporary genetic diversity and structure, indicating effects of other variable life history and ecological traits. For example, within population diversity was significantly lower and global and pairwise population measures of genetic differentiation were greater for A. spondylophylla than for A. hilliana. These differences are consistent with population genetic theory (Allendorf et al., 2013;Ellstrand & Ellam, 1993;Ghazoul, 2005;Leimu et al., 2010;Petit & Hampe, 2006) and our prediction of reduced gene flow among populations of A.
spondylophylla, which has a greater degree of geographic population disjunction, lower localized population density, smaller population size, shorter individual longevity, and lower fecundity, than A.

hilliana.
Genetic differentiation among populations located on the Hamersley and Chichester Ranges was further pronounced for A. spondylophylla, where range populations are less continuous, less dense, and smaller than those of A. hilliana (Maslin et al., 2010).  et al., 2016). Unexpectedly high measures of nuclear genetic differentiation indicate strong barriers to pollen as well as seed dispersal on more microevolutionary spatial temporal scales. Genetic structure of the nuclear genome is known to be most strongly influenced by the mating system (Duminil et al., 2007) and mechanisms of pollen and seed dispersal. Whilst Acacia typically have predominantly outcrossed mating systems (Broadhurst et al., 2008;George et al., 2008;Millar et al., 2008Millar et al., , 2019Ng et al., 2009), substantial levels of selfing have also been recorded (Coates et al., 2006;Mandal et al., 1994). Our results may be explained by increased levels of selfing or bi-parental inbreeding (mating among related individuals) although reported F IS values do not differ greatly between the study species and other Pilbara Acacia. It is possible that A. hilliana and A. spondylophylla could have more specialized insect pollinators with more limited dispersal abilities, or there could be unrecognized barriers between populations that limit the movement of important pollinators. Alternatively, it may be that stature of individual plants plays a more important role than currently recognized in long distance or localized pollen dispersal, in addition to seed dispersal, for generalist insect pollinated species in arid landscapes.
Pollen: seed dispersal ratios for A. hilliana and A. spondylophylla are most similar to those observed for the taller shrub species A. atkinsiana (12.5:1, unpublished data but see Levy et al., 2016).
However, the relative contributions of pollen and seed dispersal to overall genetic connectivity in Acacia in the Pilbara vary considerably even for tall shrub and tree species. Much higher values (43.5:1 for the tall shrub A. ancistrocarpa, unpublished data but see Levy et al., 2016), as well as much lower values (3:1 for the tall tree A. pruinocarpa, Nistelberger et al., 2020) have also been observed.
Variability in pollen: seed dispersal ratios suggest that taxon-specific life history and ecological traits other than mating system, pollen and seed dispersal syndromes, and stature also play key roles in driving diversity, differentiation, and demography. Acacia ancistrocarpa occurs in rather continuous populations at moderate density, which would theoretically facilitate greater levels of pollen dispersal (Aguilar et al., 2006;Ghazoul, 2005;Leimu et al., 2006), while A. pruinocarpa populations are more disjunct and individual trees within populations tend to be more widely spaced with possibly more limited opportunity for connectivity via pollen dispersal. More direct studies are required to confirm mating system and methods of pollen and seed dispersal in Pilbara Acacia to be able to explain this apparent variation in the relative contributions of pollen and seed dispersal between species.

| Landscape features and genetic structure
Ranges have been hypothesized to have acted as historical refugia in arid landscapes  and there was evidence of the elevated topographically complex Hamersley and Chichester Ranges acting as refugia in Eucalyptus leucophloia subsp. leucophloia . However, there was little support for this hypothesis in A. hilliana and A. spondylophylla, with no evidence for increased haplotype diversity on the Hamersley and Chichester Ranges, nor reduced diversity or signal of expansion off the ranges.
Instead, low to moderate diversity and shallow divergence among extant haplotypes, along with simple star-like networks and phylogeographic structure, support a hypothesis of relatively stable demographic history and localized population persistence both on and off the ranges. This result is concordant with scenarios of historical isolation and persistence for extant populations of other Acacia species within the Pilbara (Levy et al., 2016;Nistelberger et al., 2020). Simple star-like networks have been associated with historical population persistence and a lack of demographic signal for A. atkinsiana (Levy et al., 2016), and, while haplotype diversity and network complexity are greater, patterns also indicate historical population isolation and persistence for A. ancistrocarpa (Levy et al., 2016) and A. pruinocarpa (Nistelberger et al., 2020). A notable exception is a shared haplotype (H1) of A. spondylophylla that occurs in single populations at the north-eastern (BAT) and north-western (PAN) geographical extremes. This may be a result of homoplasy rather than long distance seed dispersal, and a similar incidence of a geographically disjunct haplotype has been noted in A. ancistrocarpa, which has a similar population range within the Pilbara (Levy et al., 2016).
Although there is general support for localized population persistence of Pilbara Acacia, greater haplotype diversity and lower phylogenetic divergence, along with a significant test of nonneutrality in extant populations, provides some evidence for historical range expansion for A. hilliana. Acacia hilliana does have greater individual longevity (50+ years vs 5-15 years) and fecundity than A.
spondylophylla (Maslin et al., 2010). Greater longevity may provide greater potential for expansion by providing greater opportunity for spatial and temporal seed dispersal. A more complex network of dispersal and significant test of non-neutrality, also suggesting a scenario of population expansion have been observed for the tall tree species, A. pruinocarpa, in the same landscape (Nistelberger et al., 2020). Acacia pruinocarpa is also long lived (50+ years) and sparsely scattered (Maslin et al., 2010), and these traits would also provide increased opportunity for spatial and temporal dispersal of pollen and seed.
Although the ranges do not appear to have had a role as historical refugia there is evidence for different patterns of contemporary genetic structure associated with these features in A. hilliana and A.
spondylophylla. Primary genetic clusters located on the Hamersley and Chichester Ranges (A. hilliana) or centrally located populations also largely on these ranges (A. spondylophylla), were differentiated from more peripheral populations at lower elevation. Limited admixture and high measures of differentiation indicate limited contemporary gene flow between these differing landscapes. Significantly greater contemporary diversity was also observed in populations located on the ranges compared to populations located in the rest of the Pilbara. A similar pattern was found for the near Pilbara endemic A. atkinsiana (Levy et al., 2016), with populations on the Hamersley and Chichester Ranges showing a greater degree of genetic connectivity and being differentiated from near coastal and other peripheral populations that occur at lower elevation. However, these patterns may be more related to taxon-specific ecological traits regarding geography and demography rather than elevation or the topographical complexity of the ranges per se. Extant populations located on the Hamersley and Chichester Ranges are less geographically disjunct and of greater local density than populations located off the ranges for all three species. Thus, the landscape scale patterns of genetic structure and diversity seen in both study species are consistent with predictions of population genetic theory that associate restricted gene flow, increased genetic differentiation and reduced diversity with geographical disjunction and low population density (Allendorf et al., 2013;Ellstrand & Ellam, 1993;Godt & Hamrick, 1993). The results strongly support a recent meta-analysis that has shown degree of population disjunction to be the greatest driver of microevolutionary population genetic differentiation in the Australian flora (Broadhurst et al., 2017).

| CON CLUS ION
Our comparative, trait-based approach suggests complex contributions of taxon-specific life history and ecological traits including species range, individual stature, and the degree of geographic disjunction in shaping genetic diversity and structure at both macroevolutionary and microevolutionary spatial and temporal scales in a common landscape. Analysis in co-distributed congeneric species provides a phylogenetically controlled evaluation of the effects of life history and ecological traits that highlights the influence of stature on genetic connectivity via both pollen and seed dispersal, a factor that has not often been recognized, as well as the effect of geographic disjunction on population genetic structure. There was little evidence for topographically complex elevated ranges acting as major historical refugia for Acacia in this landscape and the influences of life history traits are evident against a background of historical isolation and persistence of extant populations of Acacia of the arid Pilbara. Integrative studies that investigate the ecology of pollen and seed dispersal could further contribute to refining hypotheses regarding the roles of stature and localized population density in maintaining demographic and genetic connectivity in key species in arid landscapes, and other highly structured landscapes globally.

AUTH O R CO NTR I B UTI O N S
Melissa Ann Millar: Data curation (supporting); formal analysis